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Abstract — Most bounds on the stability region of Aloha give 
necessary and sufficient conditions for the stability of an arrival 
rate vector under a specific contention probability (control) 
vector. But such results do not yield easy-to-check bounds on the 
overall Aloha stability region because they potentially require 
checking membership in an uncountably infinite number of sets 
parameterized by each possible control vector. In this paper we 
consider an important specific inner bound on Aloha that has 
this property of difficulty to check membership in the set. We 
provide ellipsoids (for which membership is easy-to-check) that 
we conjecture are inner and outer bounds on this set. We also 
study the set of controls that stabilize a fixed arrival rate vector; 
this set is shown to be a convex set. 



I. Introduction 

This paper addresses geometric approximations of certain 
regions related to the stability of the slotted Aloha medium 
access control (MAC) protocol for a finite number of users. 
The stability region of a protocol is the set of arrival rate 
vectors (with elements corresponding to exogenous arrival 
rates at each user's queue) such that each queue's length 
remains uniformly bounded. This stability region is unknown 
for n > 3 users. The n = 2 user case was solved by Tsybakov 
and Mikhailov 1 1 1 in 1979, and the n = 3 user case was solved 
by Szpankowski |2| in 1994 (the expressions for ti = 3 are 
rather unwieldy). There is a large body of work giving inner 
and outer bounds on the stability region for general n (some 
highlights include IS, H, 0, ^) but many of them are 
difficult to evaluate, as we now discuss. 

Let p G [0,1]" be the contention probability vector, where 
Pi is the probability that user i will transmit a packet in any 
time slot during which user i's queue is non-empty at the start 
of that time slot, independently of anything else. Let Aa{p) 
be the stability region of the Aloha protocol associated with 
a particular choice of p, where x = {xi, . . . S A^(p) 
means the control p stabilizes each queue under arrival rates 
X. The Aloha stability region is the union of Aa(p) over all 
feasible control vectors p: 

Aa= U A^(p). (1) 

pe[o,il" 

i.e., Aa is defined by the parameterized regions Aa{p)- 

A^ = {x:3pe [0,l]":xeA^(p)}. (2) 

One limitation of many of the Aloha stability bounds in the 
literature is that they do not yield easy-to-check necessary 
and sufficient conditions for stability, meaning (to the best of 



our knowledge) the literature does not contain results of the 
desired form: 



/(x) < 

< 



X e Aa 
X ^ A^, 



(3) 



for given easily computable functions f{x.),g{x). Instead, the 
literature contains results of the form: 



/(x,p) <0 
5(x,p) <0 



X e Aa(p) 
X ^ Aa{p), 



(4) 



for given easily computable functions /(x, p), g(x, p) that 
each take as argument both x and p. These results, although 
very useful for testing stability membership in Aa(p) for a 
specific control p, are of less use in testing stability member- 
ship overall in Aa in ([T]i- Namely such results yield bounds 
of the form A^ C A^ C Ag where 



^1 



{x:3pe[0,l]":/(x,p)<0} 
{x:Vpe [0,1]" :g(x,p) <0}. 



(5) 



The main motivation behind this paper is this: most stability 
results on Aloha found in the literature are of limited value for 
evaluating stability of a rate vector x because it is not easy 
to test whether or not x £ A^ or x G Ag since the parameter 
space p e [0,1]" defining these sets is uncountably infinite. 
Consequently, the goal of this paper is to develop inner and 
outer bounds on sets like Aj,Ag that are easy to evaluate. 
The paradigmatic example we study in this paper is the set 




The expression pi Ylj^ii^ ^ Pj) is the worst-case service rate 
for user i's queue, meaning the service rate when all users 
have non-empty queues and are therefore eligible for channel 
contention. In particular, user ?'s packet is successful in such 
a time slot if user i elects to contend (with probability pi) and 
each other user j ^ i does not contend (each with independent 
probability l—pj for non-empty queue). In fact A C A^, since 
an arrival rate that is stabilized under a worst-case service 
rate is certainly stabilized under a better service rate. Further, 
for n = 2 we have A = A^ (|T|. The set A has been 
addressed in the literature in fl], fSl, pj. Massey and Mathys 
[7 1 showed A is the capacity region of the collision channel 
without feedback. Post IH showed that the complement of A 



2 



in R" is convex (thus A itself is non-convex) and characterized 
its tangent hyperplane at each point on its boundary, and 
Anantharam [4 1 showed A is the stabiHty region for Aloha for 
a certain correlated arrival process. Furthermore, A is widely 
conjectured to coincide with the Aloha stability region A a 
for general arrival processes [H, O, ||9|- Recently, by using 
mean field analysis and assuming each queue's evolution is 
independent, Bordenave et al. ifTOll were able to show that 
asymptotically this conjecture is true. 

Our approach is to provide simple geometric constructions 
that both inner bound and outer bound A, where it is straight- 
forward to check membership in both bounds (note it is not 
easy to test membership in A). Our figure of merit in evaluating 
these constructions is the volume of the set relative to the 
volume of A. In order to do so we first give in fjlll ^ new 
result characterizing vol(A). Next, in i HIIl we give a simple 
inner bound and a simple (conjectured) outer bound on A 
for arbitrary n that are inspired by the geometry of A for 
n = 2, but these bounds are loose in that their volumes 
poorly approximate the volume of A as n increases. Then jjlVI 
contains our primary result on geometric approximations of A. 
Namely, we give (conjectured) ellipsoid constructions £i,£o 
such that A/ C A C Aq where Aj ^ S\£i and Aq = S\£o, 
and iS = {x > : (x, 1) = X^iLi -^^ — ^} '■^^ simplex. 

In §V] we change our focus to a different set relevant to 
Aloha. Namely, given an arrival rate vector x, define the set of 
stabilizing controls 7'(x) assuming a worst-case service rate: 



P(x) = <^ p e [0, 1]' 



<P^Yl{l-P,), *e N >■ (7) 



This set is related to A in that 7'(x) = iff x ^ A. Whereas A 
is an important bound for the stability region A^, the set ^(x) 
is important from a more practical operational perspective: 
given a desired arrival rate vector x, it is an inner bound on 
the control options available to the network administrator. We 
establish that V (x) is a convex set. We offer a brief conclusion 
in iVi] 

II. Volume of the set A 

In this section we give an expression for the volume of A. 
Proposition 1: The set A defined in ([T]i has volume 

vol(A) = 



ke/c ^ " ' scM \n+l + T,scin]\S\ks)l 

where ("j^J^) = (fcs"sc[n]) ^ multinomial coefficient, and 

/C = i k = (fcs, S C [n]) > : |k| = ^ fcs = n - 2 

[ SC[n] 

(9) 

Proof sketch 1: Define the simplex 5 = {p > : (p, 1) < 
1} and p(p) = n"=i(l ~Pi)- One can show that the mapping 
p i-> X given by 



(p) ^PiYli^-Pj), i e W, 



(10) 



is a bijection from S to A. In fact this function is also a 
bijection from dS = {p > : (p, 1) 1} (which forms a 
"face" on the boundaiy of S) to 

dA={^:3YiedS:x,= p, - p,), ie[n]\ (11) 

o^i J 

(which forms a "face" on the boundary of A). Let J(p) = 
J(p)//9(p), where J(p) is the Jacobian for this mapping. 
The fact that dct(ayl) = a"det74 for any scalar a and any 
n X n matrix A yields det J(p) = p(p)" det J(p). Abramson 
fTTI showed that p(p)^dctJ(p) = 1 — (p, 1), which gives 
detJ(p) = p(p)""2 (1 - (p, 1)). Substituting this into the 
general expression for volume vol(A) — det J(p)dp shows 

voi(A) = n,(i - p^r-'a - E.pOdp. 

To get a better closed-form expression, we leverage results 
in IIT2I (2.3) on integration of certain functions over a simplex: 



dp = 



nn 
^=0 



(12) 



where S is the standard simplex, p = (pi,...,p„), a = 
(ai, . . . , an), and p" = I17=iPT- Defining ps = HiesP* 
the multinomial theorem gives: 



)ri-2 
-E(",')n((-i) 
keIC ^ ^ 5C[n] 



'PS 



ks 



(13) 



Next: 



n Ps" ^ nPr'"'\ «.(k) = J2ks,te [n], (14) 



SC[nl 



Substituting all these into vol(A) yields the proposition. ■ 
Note that |/C| = 1 for n ~ 2 and for any n > 2 we have 



n-2-j 



where j is the number of non-zero components of a given 
k e /C. Unfortunately \JC\ grows super-exponentially in n, 
posing huge computational overhead for even moderate n. 

III. Preliminary inner and outer bounds on A 
The Aloha stability region for n = 2 from 1 1 1 is: 



A. 



As 



{x > : 1}, 



(16) 



where the notation A^is will be explained shortly. It is straight- 
forward to establish that A^ = A for n = 2. In particular, 
for n = 2, testing membership in A = A^ is trivial since 
there is no parameterized set that must be potentially checked 
for each p G [0,1]^. Further, the shape of A suggests it 
is well-approximated by the complement of a ball, namely 
Abe = 5 \ 6(1,1), where 6(xo,r) = {x : ||x - Xo|| < r} 
is the ball centered at Xo of radius r. Fig. [T] (left) shows 



A. 



A = Aa C Abe for n = 2. 
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Fig. 1. The Asrs inner bound and A^c outer bound for n = 2 (left) and n = 
3 (right). Note the labels indicate the complement of Abe, ''(li \/n — l)n 
[0, 1]". The Asrs bound is tight for n = 2. 



These results for n ~ 2 motivate the following quantities, 
defined for arbitrary n: 



X > : ^ yi7 < 1 




{xe [0,1]": ||x-ll| >\Arn"}. (17) 

The labels "srs" and "be" stand for square root sum and hall 
complement, respectively. The latter name is appropriate since 
Abe = [0, 1]" \ &(1, V^T^), see Fig. [U 

Proposition 2: The set Ag^s in ( fTTj i is an inner bound for 
A, i.e., Asrs ^ A for n > 2. 

Proof sketch 2: Fix a point x' e Asrs- Due to coordinate 
convexity of A, Agrs, it suffices to produce a point x G A with 
Xi > x[ for each i e [n]. Set p with pi — for each 

i and set x with xi = PiWj^ii^ — Pj) for each i. Clearly 
x S A. It remains to show Xi > x'^ for each i. Note x' G Agrs 
ensures J^iPi — 1- Define independent events Ai, . . . , An 
with F{Ai) = Pi for each i E [n]. It follows that 

Then for any i, by the union bound: 

x,=p. \1-vI\JaA I >pJ 1-^P(A,)| 
Vi5^» / / V i^i / 

(19) 

and thus Xi > x[. ■ 

Conjecture 1: The set Abe in ( flTl l is an outer bound for A, 
i.e., A C Abe for n>2. 

At this point we only provide a possible methodology 
towards proving this conjecture. For any x G A we must 
show that X G Abe, meaning ||x — 1|| > -y/n — 1. To show 
this it suffices to show mirixeA ||1 — x|| > a/ti — 1, which is 
equivalent to showing 

min /(p) = ^ (l-p,n(l-?'j) ^"-1- (20) 

The approach is based on how many distinct non-zero values 
the components of p assume, i) If p has only one distinct non- 
zero value then it is easy to see /(e^) = n — 1 for each of the 



n unit vectors ei , . . . , e„ (where has a one in position i and 
zero in all the other positions). Lemma [T] states this value is 
optimal over the class of all "quasi-uniform" vectors addressed 
by the lemma, ii) If p has 2 or more distinct non-zero values 
in all its component positions, it is hoped to show that such a p 
can not be a global minimizer either because /(p) > rt — 1 or 
because p violates the Karush-Kuhn-Tucker (KKT) conditions 
required for optimality (In this case "regularity" is guaranteed, 
which justifies using KKT). 

Lemma 1: Suppose p has k non-zero elements, each equal 
to 1/fc and n — k zero elements. Then ./(p) > n — 1, with 
equality only when k ^ 1. 

Proof sketch 3: Note the aforementioned bijection property 
allows us to let J^iPi = 1- If = 1, /(p) > n — 1 is 
trivially seen to hold. If fc > 2, showing /(p) > n — 1 is 

equivalent to showing (l — A;^^)'^ ^<fc^l — \/^ ^ ij ■ 
derivative of its LHS can be shown to be negative using the 
inequality logy < y— 1. Thus the sequence {(1 — 1/fc)'^^^} is 
upper bounded by (l — i) = i. On the other hand, using 
AM-GM inequality ^1 - l/k < (1 - 1/fc + l)/2 (equality 
cannot hold when k > 2), one can see the above RHS is lower 
bounded by i. ■ 

As some further remark about the Abe bound conjecture, we 
note Cauchy-Schwarz gets close but is insufficient by itself: 



J2{i-x.)'> 



(E.(l-^.))' _ (n-E^x.)' ^ (n-1) 



The next proposition gives the volume of Asrs- Abe extends 
beyond the positive orthant so calculating its volume requires 
a difficult accounting of intersecting polar caps. 

Proposition 3: The volume of the inner bound Agrs is 



vol(As 



2" 



< vol(A). 



(21) 



Proof sketch 4: Use change of variable yi = y/x%i for i G 
[n], and define Sk = J^i^k Hi ^^^^ '■^^ volume becomes 

vol(Asrs) = 2" / ls^<iyidyi ■ ■ ■ yndyn- (22) 

■/[O.l]" 

Integration gives: 

/o "^'""^((1 - ^fc+i) - Vkr-'vkdyk i-iy f2k - 2 



(l-S,+ i)2'= 



E 



i + 2 V i 



(23) 

for each fc. By induction on ( |22] ) and ( |23] ) we can estabUsh: 
v„,(A..,,^2"n(gti|(^\-^)). (24) 



It is straightforward to show that 

2fc-2 ^_^y /2fc„2 



E 



- * + 2 V ^ J (2fc)(2fc-l)' 11 



1 

Y[{2k){2k^l) = {2n)l 
(25) 



As shown in Fig. [3] vol(Asrs) scales rather poorly in n 
relative to vol(A), therefore we seek improved bounds for A. 
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IV. Ellipsoid inner and outer bounds on A 

The fact that the Abe curve is superior to the Agis curve 
in Fig. [3] suggests the complement of A may be well approx- 
imated by appropriately chosen ellipsoids. To this end, we 
consider ellipsoids £ of the form (ifTsl): 

(26) 



£ = {x : (x - c)"rR-i(x - c) < 1} 



Here c is the center of the ellipsoid and the n x n matrix 
R has the spectral decomposition R = QDQ^^ where Q = 
[qi • • ■ Qn] is orthonormal and holds the eigenvectors of R 
which are the directions of the n axes of the ellipsoid, and 
D = diag((ii, . . . ,dn) holds the eigenvalues of R, and each 
o-i = \fdi is the semi-axis length in the direction q^. 

We assert that the symmetries in the set A suggest sym- 
metries for £. In particular, the center c should lie on the 
ray (1,...,1), i.e., c = cl = (c, ...,c) for some scalar 
c > 1/n to be determined. Further, one axis of the ellipsoid 
should point along this same ray, i.e., have associated unit 
vector qi = 1 = (1, . . . , l)/-yn and associated length a\. 
Imagine we rotate the coordinate system so that the vector ei 
is aligned with the vector 1. Then this induces a permutation 
symmetry among the remaining rotated n — 1 coordinates. 
Hence, the ellipsoid is spherical in those directions (meaning 
a2 = • • • = a,i), leaving three degrees of freedom: c, ai, 02- 

Our approach is to approximate the surface 9A with part of 
the surface of an ellipsoid placed in the complemenet of A, and 
then form inner and outer bounds on A by subtracting these 
ellipsoids from the unit simplex S. The following conjecture 
defines the proposed ellipsoids and the corresponding regions 
that we believe to be inner and outer bounds on A. 

Define the ellipsoids £/ , Eq with center c = (c, . . . , c) and 
R = QDQ^^ with orthonormal Q — (qi, . . . , q„) for qi = 
(1, . . . , l)/^/n and q2, . . . , q„ otherwise arbitrary. The two 
distinct semi-axes lengths are: 

2 



£1 ■■ 

£0 



1 



1 \ ^ 1 



af = {nc — l)c, 



^2 _ 
2 na^ — (nc— 1)^ 

flj = (ri — l)c 



(27) 

Define also the regions Aj = S \ £1 and Aq = S \ £o- 

Conjecture 2: The regions A/,Ao form inner and outer 
bounds on A: A/ C A C Aq for n > 2. 

The expressions for ai, 02 for £j, £0 are derived essentially 
based on Lemma|2] Corollary [T] and Lemma|3]given below. In 
words, £j is such that it passes through ei, . . . , e„ and the "all- 
rates-equal" point m = ml e dA where m — i(l — ^)"^^- 
It can be shown that d£i is tangent at m with dA. £0 is such 
that it passes through each and further d£o is tangent with 
dA at each e^. Note when c = 1, Aq = Ai,c since in this case 
£0 is a ball with ai = 02 = %/ n — 1 . 

The following lemma says that, under given conditions, the 
ellipsoid is invariant to the choice of q2, . . . , qn provided Q 
is orthonormal and qi = l/\/n. 

Lemma 2: For any ellipsoid in the form of ( |26] |. if qi = 
{\, . . . ,1)1 \/n and D = diag(af, a|, . . . , Oj), then R^^ = 



Cl + 03 I, where ( = ^ (oi ' 



1 is the n X n matrix 



with lij = 1 for each i,j and I is the identity matrix. 

The proof is straightforward and omitted. This simple form 
of R~^ allows us to characterize the ellipsoid that passes 
through each e^: 



Corollary 1: For any c > 1/n and ai > y/n{c — ^/n), 
setting 

2 _ {n-l)al 

- — 2 7 TTI (28) 

rta| — [nc — 1 

ensures G d£ for i G [n]. 

The next lemma tells what is needed for d£ and dA to share 
a common tangent point: 

Lemma 3: Define qt = 1 — pt, d£ and dA share a point 
of tangency at xt if for each i G [n — 1]: 

'a2\^ _ ilt.i - qt,n)^^=iXt,j +n{qt,nXt,t - qt.iXt.n) 

/ {qt.i - qt.n){^]=iXt,j - nc) 

(29) 

Proof sketch 5: Using the bijection ( fTOl ) between dA and 
dS it follows that I?] there is a unique pt €E dS with (p, 1) — 
1 corresponding to x^. Post [H established that the tangent 
hyperplane to A at a point Xt G i9A with Xt 1— > pt is (1 — 
Pt,x) = p(pt), for p{pt) = 11^(1 - Pt,i)- Next, using the 
implicit function theorem it is straightforward to establish that 
the tangent hyperplane to an arbitrary ellipsoid £ ( |26l ) at a 
point xj e d£ is given by (w,Xt) =6 where 

R-i(xf - c) 



-, = (w,Xi) 



(30) 



R„ ^ (xf - c) 

Here R~^ is the n*'' row of R^^. Finally equating (after 
appropriate normalization) the components of the normal 
vector of dA with those of d£ yields this lemma. ■ 

Fig. |2] illustrates the ellipsoids £i,£o and the associated 
inner and outer bounds A/, Aq for the special case of n = 2 
for c = 0.55 and c = 2. It is observed that the accuracy of the 
bounds improves as c is increased. 
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Fig. 2. The (conjectured) inner (left) and outer (right) bounds on A using 
c = 0.55 (top) and c = 2 (bottom). The quality of the bound improves as c 
increases. 

If the conjectured bounds hold then dA is increasingly 
tightly "sandwiched" between part of d£i and d£o as c — )■ oo. 
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This sandwiching is asymptotically tight at the all-rates- 
equal point m (and by construction always tight at each 
ei). Although this tightness decreases in n it is bounded. 
Specifically, the ellipsoids £i,£o viewed in the limit as 
c — )• oo have axes ratios given by limc^oo 1/0,1 



^ Q = 1 and 



limr 



for each n (see (|27])). 



^2.0 - 2(i_(i_i)"-i) 

This latter expression equals 1 when n = 2 and monotonically 
decreases in n to e/(2{e— 1)) « 0.791. Thus when n — 2, the 
£i,£o are asymptotically equal as c — > 00. For large c and for 
arbitrary n, the axis ratio for ai is always approximately one, 
while the axis ratio for is bounded between (0.791, 1.0). 




Fig. 3. The volumes of A, Asrs, Abe, A/ , Aq versus n. 

Fig. [3] shows the volumes of A, Agis, Abe, A/, Ao versus 
n = 2,..., 6, with c — 100. Monte-Carlo simulation was 
used to obtain the volumes of Abe, A/, Aq, while ([8]) and ( ITTT i 
were used for the volumes of A and Asis, respectively. 

V. Stabilizing controls for a rate vector 

We now change gears and discuss the set 7'(x) in (|7|, which 
is the set of contention probabilities that can stabilize the given 
arrival rate vector x assuming worst-case service rates. The 
motivation is that knowing a rate vector is stabilizable (i.e., 
whether or not x e A or x G A^) is useful only if you 
also know how it may be stabilized, i.e., 'P(x). Note 7^(x) C 
Va (x) for Va (x) the set of controls that stabilize Aloha under 
arrival rates x. Our next result establishes that 7'(x) is convex. 
As an example, the region 7^(x) for x = (0.192,0.052) is 
shown in Fig. |4] 

Proposition 4: The set 7'(x) in (|7]i is convex for all x G A. 

Proof sketch 6: We can decompose 7'(x) as the intersection 
of the n regions V{x) — PlILi 'Pii^i) where 



i € n 



(31) 



'Pi{xi) = <P ■■ Xi < p, J|(l - Pj] 
[ J^'. 
Since convexity is preserved under intersection, it suffices to 
show that each Vi{xi) is convex. Letting p_i be the (n — 1)- 
vector p\pi and /i(p-i) = Xi/ Oj^ill ~-Pj)' we can express 

V^{xi) = {{v-uPi) ■■ < pJ , « e [n]. (32) 

But ( |32] | is the epigraph for /; and thus Vi{xi) is a convex 
set iff fi is a convex function (|13| p. 75). The function 

/i(P-i) = log/»(p-i) = logxi + ^ - log(l -pj) (33) 



is a convex function of p_,: since — log(l — Pj) = — log(l — 
ejp_i) is a convex function of p i and fi is therefore a non- 
negative weighted sum of convex functions, and is therefore 
itself convex ( ifTSll p.79). Since fi is log-convex it is also 
convex (fT3l p. 104). ■ 




Fig. 4. The region P(x) for x = (0.192, 0.052). 

VI. Conclusion 

We have provided easy-to-check inner and outer bounds 
on the set A, which is a paradigmatic example of an inner 
bound on the Aloha stability region that is difficult to check 
membership. There are many directions to pursue to extend 
this work, most notable is to prove the Conjecture |2] (which 
extends Conjecture [U. Explicit expressions for the volumes of 
these sets as a function of n are desirable, but appear difficult. 
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